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ABSTRACT 

The 'Internal Linear Combination' (ILC) component separation method has been ex- 
tensively used on the data of the WMAP space mission, to extract a single component, 
the CMB, from the WMAP multifrequency data. We extend the ILC approach for re- 
constructing millimeter astrophysical emissions beyond the CMB alone. In particular, 
we construct a Constrained ILC to extract clean maps of both the CMB or the ther- 
mal Sunyaev Zcldovich (SZ) effect, with vanishing contamination from the other. The 
performance of the Constrained ILC is tested on simulations of Planck mission obser- 
vations, for which we successfully reconstruct independent estimates of the CMB and 
of the thermal SZ. 

Key words: Cosmic Background Radiation - Methods: data analysis 
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1 INTRODUCTION 

The separation of components in observations of the Cos- 
mic Microwave Background (CMB) is an important part of 
the processing and analysis of such observational data. Var- 
ious component separation methods have been developed to 
extract the emission of a single component (or of several 
of them) out of multifrequency observations (see, e.g., De- 
|labrouille fc Ca rdoso (2009) for a review). 

Often, such methods assume that the observations are 
a linear mixture of unknown components (or sources), in 
which case the data are modelled as: 



Xiip) 



(1) 



where Xi(p) axe N ^ a observed maps (p indexing the pixel), 
Sj(p) are iV comp templates for unknown components of in- 
terest, and rii(p) are maps of noise for the observations. The 
coefficients Aij define a mixing matrix. 



after it has been determined using one of the forementioned 
blind component separation methods), the actual compo- 
nent separation is solved by inversion of the linear system of 
equation [T] Methods for doing so in the presence of instru- 
mental noise have been investigated by a number of authors 



ITegmark fc Efstathiou 1996 Bouchet fc Gispert 1999 



Hobson et al. 1998 Delabrouille et al. 2002 1 



( Delabrouille et al. 



Blind component separation methods such as SMICA 

[20031 iCardoso et al.l [2008]), FastICA 1 2004 1 |Park et al.||2007||kim et al.||2008||Delabrouille et al 



In CMB observations, in practice, at least one of the 
columns of the mixing matrix is usually known (i.e. one 
mixing vector) - that of the CMB. On the other hand, some 
components cannot be modelled as a single template which 
is simply scaled by mixing coefficients (e.g. emissions from 
the galactic interstellar medium). Then, the data model of 
equation [I] fails, and one has to resort to other approaches. 
For these reasons, the so-called 'Internal Linear Combina- 
tion' or ILC, wich does not assume any particular parametri- 
sation for foreground emisson, has been extensively used in 
the analysis of the maps obtained by the WMAP satellite 
to extract a CMB map (jBennett et al. 2003[ |Eriksen et"aT 



( Hyvarinen[ |1999| |Maino et al.| |2002| )7JA DE (|Cardoso 



1998), CCA (Bonaldi et al. 20061 or CMC A (Bobin et al 



2009) 



2008) are primarily designed to solve the problem of esti- 



mating the sources S;, separated from one another, in the 
case where the observations can be modelled as in equation 
[I] with the matrix A unknown. In practice, the first (and 
most difficult) task is to estimate A. 

Then, if/when matrix A is known (either a priori, or 
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The Planck mission, launched on the 14th of May 2009, 
is a third generation CMB experiment. It observes the sky 
in 9 frequency channels ranging from 30 to 857 GHz. The 
high frequency instrument of Planck ( jLamarre et al.[ |2000[ 
20031, in particular, has been designed with bands centred 



at the minimum, the zero, and the maximum of the thermal 
SZ emission. The extraction of clean CMB and SZ maps, 
of a catalog of galaxy clusters selected by their thermal SZ 
effect, and the investigation of bulk flows in the large scale 
velocity field from kinetic SZ effect towards galaxy clusters, 
are part of the scientific programme of Planck. 
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For these projects, accurate separation of CMB (and which can be recast as 



kinetic SZ) from thermal SZ is important (Aghanim et al 



2001 


2003). 


CMB 


maps 



2001 20031. In particular, residuals from thermal SZ in the 



They can also bias the estimation of cosmological param- 
eters, in particular those which depend drastically on the 
small scale CMB power spectrum, as the scalar spectral in- 



2009) 



dex and the reionisation optical depth ( Taburet et al. 

In the present paper, we address the problem of ex- 
tending the ILC method to separate several components of 
interest with known "mixing vector" in multifrequency ob- 
servations such as those of Planck. The method, denoted as 
Constrained ILC, is of interest for separation of CMB and 
thermal SZ components with vanising residual contamina- 
tion of one by the other. 



2 ILC ESTIMATION OF CMB AND SZ 
2.1 Astrophysical components 

Astrophysical components can be separated into two broad 
categories: diffuse components and point sources. 

Point sources (i.e. unresolved objects) are typically de- 
tected and identified with specific methods, based on spatial 



filtering (see, e.g., Barreiro ( |2009[ ) for a review). Such meth- 
ods are effective when detectable sources are isolated. Ex- 
periments observe also a background of faint sources which, 
due to limited sensitivity and resolution, cannot be detected 
individually. Such a background is treated as a diffuse com- 
ponent. 

Other diffuse components of interest comprise the CMB 
and diffuse emission from the ISM. SZ effects (thermal and 
kinetic), can be considered either as point sources or as dif- 
fuse components, depending on the resolution, sensitivity, 
and frequency coverage of the experiment considered. With 
Planck, a number of observable clusters will actually be re- 
solved, and we will consider here the thermal SZ as a diffuse 
component. 

The emission of the Interstellar Medium (ISM) is com- 
plex, as it involves several processes of emission which are 
not fully independent. Hot gas emits through synchrotron 
radiation, warm gas by free-free emission, cold dust by grey- 
body emission, and plausibly by dipole emission from ro- 
tating dust grains. Molecules emit through molecular tran- 
sitions. The dependence of these emissions with frequency 
is a function of additional parameters which vary over the 
regions of emission, e.g. electron temperature for free-free, 
distribution of electrons as a function of energy for the syn- 
chrotron, composition and temperature of molecular clouds 
for molecular line spectra... In addition, as all the ISM is con- 
centrated towards the galactic plane, all of these processes 
do not result in independent, nor even simply uncorrelated, 
emissions. 



2.2 The ILC method 

The standard ILC assumes very little about the model of 
the data. It simply assumes that all available maps (N b B 
maps, indexed by i) can be written, for all pixels p of the 
observed maps, as 



x(p) = as(p)+n(p), (3) 

where x(p) is the vector of observations (N h s maps), s(p) a 
single map of a component of interest, a is a known mixing 
vector which does not depend on p, with as many entries as 
there are channels of observation, and n includes instrument 
noise as well as all other astrophysical emissions. For all 
channels (iV bs maps), it is assumed that all observations are 
at the same resolution, although a harmonic or needlet space 
implementation of the ILC permits to deal with channel- 
dependent resolution in a simple way (assuming symmetric 
beams). 

The ILC provides an estimator silc of s by forming 
the linear combination s(p) = w t x(p) of the observed maps 
which has unit response to the component of interest (i.e. 
w t a = 1) and has minimum variance. Straightforward al- 
gebra leads to ILC coefficients such that silc of s is given 
by: 

e (4) 



silc = 



a' R- 



where R is the empirical covariance matrix of the observa- 
tions. 

The ILC component separation method has advantages 
and drawbacks. The main advantage is that it does not as- 
sume a model for the components we are not interested in, 
i.e. all the components whose contributions are collected into 
a single nuisance term n(p). Drawbacks include the existence 
of a bias induced by any empirical correlation between the 
component of interest and the nuisance term, as described 



in the appendix of Delabrouille et al. ([2009]), and the need 
to know the coefficients a* with some accuracy, especially 



for sensitive experiments (Dick et al. 2010) 



Note that the ILC relies on the component of interest to 
be uncorrelated with the contaminants, i.e. (s(p)ni(p)) — 
for all channels of observation i. 

In its simplest implementation, the ILC is performed on 
the complete maps, and one single global matrix R is used. It 
is possible, however, to decompose the original maps as sums 
of different data subsets, covering each a different region 
in pixel space or in harmonic space, to apply independent 
versions of the ILC to the different data subsets, and then to 
recompose a map from all these independent results. In the 
present paper, all the experiments done on simulated maps 
are performed in harmonic space. We have R = H(£), and 
each such covariance matrix is estimated as 
i+Ai 



RyW 



1 

iV> 



E E 



(5) 



where the spherical harmonic coefficients of map 

Xi(p) and where Ne is the number of modes (£,m) in the 
window [£- A£,£ + M\. 



N e = (£ + A£ + if -(£- A£)' 



(6) 



Xi{p) = CLiS(p) + Ui(p), 



(2) 



Each single K(£) is obtained as the average over a window 
in t. 

2.3 ILC and the SZ effect 

So far, the ILC has been used almost exclusively to extract 
a CMB map, but it can be used in a similar way to extract 
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any single component which is described as a single tem- 
plate scaling with frequency, provided that the appropriate 
column of the mixing matrix is known, and that this tem- 
plate is not correlated with other emissions present in the 
data set. 



Sunyaev-Zel'dovich emission (Sunyaev & Zeldovich 



19721, which can be observed towards clusters of galaxies 



at millimeter wavelengths, arises through inverse Compton 
scattering of CMB photons off hot electrons of the intra 
cluster gas (see, e.g. Birkinshaw ( 1999 1 or Carlstrom et al 
(20021) for a review). 



The thermal SZ emission, in the non-relativistic approx- 
imation, is given by the product of a template map yip) (the 
map of the cluster Compton parameter, proportional to the 
integral over the line of sight of n e T e ), and of a known emis- 
sion law f(f). The Compton parameter is given by: 

V = I %rn e a T dl, (7) 

il.o.s. "»eC 2 

where k is the Boltzmann constant, m e the electron mass, 
c the speed of light, cry the Thomson cross section, and n e 
and T e the electron density and temperature respectively. 
The frequency dependence of the effect is: 

B =M + i 



with 



/(f) =*(")- 



,(*) = 



hv 



(8) 



kTcus 

The coefficients <n of equation[2] for the thermal SZ, are the 
integral of the emission law f(u) in the frequency bands of 
the instrument used to observe the sky. 

Kinetic SZ emission has the same emission law as CMB 
anisotropics, and is proportional to the peculiar velocity of 
the scattering electron gas (i.e. of the cluster) along the line 
of sight. Although the kinetic and thermal SZ effects arise 
from the same set of galaxy clusters, their correlation van- 
ishes because of the sign dependence of the kinetic SZ tem- 
perature. 

Hence, it is possible to extract a map of the thermal 
Sunyaev Zel'dovich (SZ) effect from multifrequency observa- 
tions with an ILC in the same way as is done for the CMB. 
This method has been implemented and tested in data chal- 
lenges using simulated data sets, organised in the context of 
the preparation of the Planck mission ( |Leach et al 



Melin et al. 2010 1 



2008 



Note that on the other hand, a CMB map obtained by 
an ILC contains emission from both primary CMB anisot- 
ropics and kinetic SZ effect. The latter is typically very small 
compared to primary CMB anisotropics. For high resolution 
high sensitivity maps however, the kinetic SZ can be sepa- 
rated from the CMB by matched filtering, using the cluster 
profile estimated from the thermal SZ map. 

2.4 A two-component model for the ILC 

Without assuming much about the detailed properties of the 
emissions of other foregrounds, and supposing we are inter- 
ested mostly in the CMB and the SZ, both thermal (tSZ) 
and kinetic (kSZ), the observational data can be modeled 
by an extension of equations [2] and [3] as: 

xi(p) = <Hs(p) + hy(p) +m(p) (9) 



zip) = as(p) + by(p) + n(p) 



(10) 



where, as in equation [3J s(p) is the CMB map (includ- 
ing kSZ), but where now y(p) is the thermal SZ map and 
n = {n.i(p)} includes both instrumental noise and unmod- 
eled astrophysical foregrounds (i.e. all sky components ex- 
cept CMB, kSZ and tSZ), in all frequency channels. Vector 
a = (1, 1, 1)* is the CMB mixing vector and b is the mix- 
ing vector of the thermal SZ, as derived from Eq. Q. 



2.5 Independent ILC estimation of CMB and SZ 

It is a straightforward process to reconstruct both a 
(CMB+kSZ) and a tSZ map independently by implement- 
ing two separate ILC, one for the (CMB+kSZ), and one for 
the tSZ. Standard ILC gives: 



s(p) 



a* R- 



x(p), 



yip) 



b f R- 



■ x{p). 



a'R-ia " ""' f/R-ifo 
Those estimates, however, do not take fully into account 
the prior knowledge about the existence of two components 
entering in the observations as described by equation [9] In 
particular, the ILC weights used for CMB reconstruction 
do not guarantee that the reconstructed CMB contains no 
thermal SZ. Similarly, the weights used to reconstruct the 
thermal SZ do not guarantee that the reconstructed thermal 
SZ does not contain any CMB (and kinetic SZ). In fact, using 
the solutions above with the data model of equation [9] we 
have: 

of R ^b of R 1 . . 

s(p) = s[p) + ^ y(p) + ^ n(p) 



a' R- 



a* R- 



The second term on the right hand side is the contamination 
of the recovered (CMB+kSZ) by thermal SZ. Similarly, when 
one recovers a tSZ map, we get: 



yip) = yip) + 



b'R-ifo " b t R~ 1 b 
These solutions minimise (by construction) the total vari- 
ance of the reconstructed maps, but both maps contain con- 
tamination from the other component of interest. For certain 
applications, this contamination is not acceptable. For in- 
stance, significant thermal SZ leaking into the reconstructed 
CMB will make it difficult to extract any kinetic SZ infor- 
mation from the reconstructed (CMB+kSZ) map. 

2.6 Constrained ILC estimation of CMB and SZ 

Now assume that we want to ensure that each reconstructed 
map of interest contains no contamination from the other 
component, i.e. the CMB map contains no contribution from 
the thermal SZ, and the thermal SZ no contribution from 
the CMB (nor from kinetic SZ). 

In the spirit of the ILC, we then look for a minimum 
variance estimate s of the CMB marQ s as a linear combi- 
nation of the frequency maps a^: 

s = w t x 

1 The same derivation can be performed for thermal SZ recon- 
struction, by simple inversion of the roles of a and b and of s 
and y. 
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0.0 0.10 mK CMB -0.060 ^^^^n 0.060 raK CMB -0.40 0.40 mK CMB 

(33.0, 8B.O} OoJnctic (33.0, 8B.0} OoJnctic (33.0, 8B.0} OoJnctic 



Figure 1. Simulated observations: A patch of the simulated sky located at high galactic latitude, around galactic coordinates of 
(7, ft) = (33°, 89°). From top left to bottom right: CMB (including kinetic SZ), Kinetic SZ, Thermal SZ, other foregrounds (at 143 GHz), 
Noise of the 143 GHz channel, Total observed emission (including noise) at 143 GHz. 



where the weights Wi have to satisfy the constraints: 

w t a = l, (11) 

w>*b = 0, (12) 

so that we conserve the CMB signal (including kSZ) and 
completely eliminate the thermal SZ signal. The weights 
thus are solutions of 



d 
dun 



\w Rw + A (l — w o) — ftw'b 







where R is the empirical covariance matrix of the observed 
maps and A, n are Lagrange multipliers. We find the solu- 
tion: 



w — AR 1 a + ^iR b. 



(13) 



Applying constraints (111, ( 12 1 to the solution ( 13 1, we have 
to solve the system 



Aa'R 1 a + fia t R 1 b = 1, 
\b t R- 1 a + [xb t R- 1 b = 0. 



(14) 
(15) 



If b is not proportional to a (which is guaranteed since the 



CMB and the thermal SZ do not have the same emission 
law), there is a unique solution to the system, for which the 
Lagrange multipliers are: 



A : 



M = 



WR^b 



a^R^b 

tR-l^ ( f>*R-l 



a*R-!a lb t R- 1 b) - a'R-ife 



2, (16) 
— • (17) 



The weights of the constrained ILC (for CMB+kSZ recon- 
struction) are given by 



b'R^b ) a'R- 1 - [ a'R^b ) b'R" 1 



a'R-ia] ( b*R _1 b 



^R-^y 



(18) 



A similar expression with the role of a and b exchanged 
holds for the weights to be used for estimating the thermal 
SZ map y(p). 

This is the generalisation of the ILC component sepa- 
ration when two components are recovered, and when one 
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Recovered CMB+Kinetic SZ: Std. ILC Reconstruction error: Std. ILC Reconstruction error: Constr. ILC 
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Figure 2. CMB extraction: Left panel: The CMB reconstructed by standard ILC from simulated Planck observations (visualisation of 
the small patch of figureJTJ. Reconstruction looks excellent. Middle panel: a detailed examination of the reconstruction error (recovered 
CMB minus input CMB, including the kSZ), reveals dark patches towards the direction of galaxy clusters. This is due to the leakage of 
thermal SZ into the reconstructed CMB. Right panel: a constrained ILC guarantees the absence of contaminating tSZ in the reconstructed 
map, with minimal impact on the total error variance. 



imposes that there is no leakage of one component in the 
reconstructed map of the other. 



2.7 Application to simulated Planck observations 

We investigate the performance of standard ILC and con- 
strained ILC for separating the CMB and the thermal SZ 
using observations such as those of the Planck mission. 



2. 7. 1 Simulations 

Our investigations are carried out on sky simulations gen- 
erated with the Planck Sky Model (PSM) version 1.6.3 
for all Planck LFI and HFI channels. Sky simulations in- 
clude Gaussian CMB generated assuming a Ct model fit- 



ting the WMAP 5 year observations ( |Hinshaw et al. 2009 ) , 
thermal and kinetic SZ effect, four components of galac- 
tic ISM emission including thermal and spinning dust, syn- 
chrotron, and free-free, and emission from point sources (ra- 
dio and infrared). The resolution and noise level of the ob- 
servations correspond to nominal mission parameters as de- 
scribed in the Planck 'Blue Book'. Some details about PSM 
simulations can be found in |Leach et al. ( 2008 1 and Be- 
|toule et al.| ( |2009| . Figure [T] displays the simulated emis- 
sion in a small patch centred around galactic coordinates 
of (I, b) = (33°, 89°) (selected, in our full sky maps, for the 
presence of a bright galaxy cluster at high galactic latitude, 
where galactic foregrounds are low, which permits to illus- 
trate best the separation of CMB and SZ). 



2.7.2 ILC results: the CMB 

As a first step, we implement the usual single component 
ILC (standard ILC) for both the CMB and the thermal SZ. 
For this particular application, a small galactic mask is used 
to blank out the region of strongest galactic emission, which 



permits to implement the ILC in harmonic space (i.e priority 
is given to the localisation of the filters in harmonic space, 
rather than pixel space). The spectral statistics K(£) are 
computed in windows of I of width Al — 50 at low £, and 
A£ — 20 for the highest multipoles. The ILC is performed 
independently for each i. 

Figure [2] displays the results of the reconstruction of the 
CMB. The left panel shows the standard ILC reconstruction 
of the CMB map, which can hardly be distinguished from 
the left panel of figure [I] (the reconstruction is visually ex- 
cellent, albeit one may notice faint small scale granularity 
caused by noise). The middle panel shows the error (differ- 
ence input-output) map using a standard harmonic domain 
ILC. Negative patches, which could not be picked out in the 
reconstructed CMB map, are clearly seen in the direction of 
the brightest clusters : the standard ILC, clearly, does not 
reject perfectly the thermal SZ effect. The amplitude of the 
thermal SZ leaking into the map is of 0.1 mK (thermody- 
namic) for the brightest cluster in the field, about 2.5 times 
the kinetic SZ effect for that particular galaxy cluster (and 
well above the Planck noise level). 

The right panel shows the reconstruction error when 
the constrained ILC is used instead (the corresponding CMB 
map cannot be distinguished by visual inspection from what 
is obtained with the standard ILC, and is not displayed. As 
expected, there is no indication in the error map of leakage of 
thermal SZ into the reconstructed map of CMB + kinetic SZ. 
In the present case, the impact of the additional constraint of 
vanishing thermal SZ contribution in the CMB has negligible 
impact on the total noise level. 

That example demonstrates that with data sets such as 
those of Planck, constraining the ILC to reject the thermal 
SZ effect makes it possible to avoid contaminating the CMB 
with SZ, with very little impact on the overall level of noise 
in the reconstructed CMB. This statement, however, does 
not hold for any possible experiment. In the case of obser- 
vations with more noise and less frequency channels, it may 
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Reemrei-ed Thermal SZ: Std. ILC 




ulations of Planck observations, and have shown that they 
permit to reconstruct both a CMB and a thermal SZ map 
with excellent performance. 
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Figure 3. Thermal SZ standard ILC reconstruction centred 
around galactic coordinates of (I, 6) = (33°, 89°) (5' resolution). 

well be that constraining the ILC would actually result in 
significantly increased noise level, for negligible gain on the 
contamination by SZ. 

2. 7.3 ILC results: the SZ effect 

Similarly, we show in figure [3] the reconstruction of thermal 
SZ effect by a standard ILC in that same area of the sky. 
Clusters are clearly visible, which confirms the sensitivity of 
Planck for detecting galaxy clusters and doing SZ science. It 
is clear that the very specific SZ emission law, with negative 
signal below 217 GHz and positive signal at higher frequen- 
cies, helps to separate it from other emissions effectively. 

We compare the performance of the Standard ILC result 
(shown in figure [3]l , and a constrained ILC reconstruction 
(not shown), where the solution is constrained to perfectly 
reject components with the same color as that of CMB an- 
isotropics (i.e. both the CMB and the kinetic SZ). The result 
of the constrained ILC is not displayed since it is visually 
indistinguishable from the standard ILC result. Visual in- 
spection of reconstruction errors (reconstructed thermal SZ 
minus true input thermal SZ at the same resolution) does 
not reveal any particular feature connected to any astro- 
physical component, in either reconstructed SZ map. 

Hence, in the case of the reconstruction of the thermal 
SZ, the constrained ILC result is very similar to the standard 
one. The main reason for this lack of qualitative difference 
between the two filters is simple: as the CMB is seen by all 
channels with very good signal to noise ratio, the standard 
ILC always adjusts the ILC weights to null-out the CMB 
(or nearly so). Little is gained by imposing this constraint 
explicitely a priori. 



3 CONCLUSION 

In this article we have developed a Constrained ILC method, 
and have shown how it can be used to reconstruct the CMB 
and the thermal SZ component with vanishing contamina- 
tion of one by the other. We have applied the filters to sim- 
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